Anomalous temperature dependence of the Casimir force for thin metal films 
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Within the framework of the Drude dispersive model, we predict an unusual non-monotonous 
temperature dependence of the Casimir force for thin metal films. For certain conditions, this 
force decreases with temperature due to the decrease of the metallic conductivity, whereas the force 
increases at high temperatures due to the increase of the thermal radiation pressure. We consider 
the attraction of a film to: either (i) a bulk ideal metal with a planar boundary, or (ii) a bulk 
metal sphere (lens). The experimental observation of the predicted non-monotonous temperature 
dependence of the Casimir force can put an end to the long-standing discussion on the role of the 
electron relaxation in the Casimir effect. 

PACS numbers: 11.10. Wx, 73. 61. At, 



The Casimir effect is one of the most interesting macro- 
scopic manifestation of the zero-point vacuum oscilla- 
tions of the quantum electromagnetic field. This effect 
manifests itself as the attractive force arising between 
two uncharged bodies placed in the vacuum due to the 
difference of the zero-point oscillation spectrum in the 
absence and in the presence of them (see, e.g., the mono- 
graphs [H, and review papers [1, 3). 

The Casimir effect attracts considerable attention be- 
cause of its numerous applications in quantum field the- 
ory, atomic physics, condensed matter physics, gravi- 
tation and cosmology 0, 0, Q, [||. The noticeable 
progress in the measurements of the Casimir force Q 
has opened the way for various potential applications in 
nanoscience , particularly, in the development of nano- 
mechanical systems 0, . 

In spite of intensive studies on the Casimir effect, it is 
surprising that such an important problem as the temper- 
ature dependence of this effect is still unclear and is still 
an issue of lively discussion 0, 0, [l(| . The central point 
in this discussion is if the Lifshitz formula (sec, c.g 
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is applicable or not for lossy media. The authors of Rcf. 8 
have argued that the Drude dispersion relation for a lossy 
medium leads to inconsistencies because the reflection co- 
efficient tte for the TE electromagnetic mode becomes 
discontinuous when the imaginary frequency £ = — iu> 
tends to zero. Therefore, instead of the Drude dispersion 
relation for the high-frequency dielectric permittivity e, 



e(*C) = 1 + 



(1) 



where uj p and v are the plasma frequency and the re- 
laxation frequency, authors of Ref. [s| suggest the same 
equation, but with v = 0. Bostrom and Sernelius [t| 



have been the first to inquire whether this prescription is 
correct. They argued that in view of a realistic disper- 
sion relation, the TE mode should not contribute to the 
Casimir force at zero temperature. Later, the authors of 
Refs. 0, fTol have shown that the mentioned discontinuity 
of tte at ( — > does not lead to any physical difficulty 
or ambiguity. 

For the case of zero temperature, the essence of the 
problem can be reduced to the following fundamental 
question: Can the Casimir force "feel" the dissipation 
parameter (the relaxation frequency v) at zero temper- 
ature when the dissipation itself is absent? According 
to Ref. la, the answer is "no" . However, the authors of 
Refs. 0,|l3 conclude that the answer should be "yes" . 

In this paper, we pay attention to an important feature 
of the Casimir force that can be demonstrated within 
the frame of the Drude dispersion model. There exist 
two competing phenomena that obviously determine the 
temperature dependence of the Casimir force. On the one 
hand, an increase in temperature leads to an increase of 
the relaxation frequency and, therefore, to a decrease of 
the metal conductivity and to a decrease of the Casimir 
force. On the other hand, when increasing the temper- 
ature, the Casimir force increases due to the growth of 
the thermal radiation pressure. The competition of these 
two effects can result in a non-monotonous temperature- 
dependence of the Casimir force. 

The experimental observation of such an anomalous 
temperature dependence of the Casimir force might be 
a direct justification of the applicability of the Drude 
model. However, this temperature effect for bulk metals 
is very difficult to observe because of its small magni- 
tude. Indeed, the relative contribution |F ra( j/F(T = 0)| 
of the thermal radiation force F ra d into the Casimir force 
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F is proportional to (kTa/hc) 4 <C 1 for realistic dis- 
tances a <C clt = hc/kT between bulk metals. Here k 
is the Boltzman constant and c is the speed of light. The 
temperature-dependent part of the term F v , related to 
the relaxation frequency is, is very small because it is 
proportional to the small surface impedance of a metal. 
Therefore, for bulk samples, the Casimir force is observed 
to be slowly increasing with T due to an increase of the 
radiation term i^ad- 

In this paper, we predict a decreasing Casimir force 
with T and show that the difficulties mentioned above, 
for the observation of the anomalous temperature depen- 
dence of the Casimir force, can be significantly dimin- 
ished if we consider the interaction of thin metal films, 
instead of only between bulk samples. As was derived in 
Ref. Il2j, the temperature effects in the Casimir force can 
be brought to the forefront if the film thickness d is the 
smallest parameter with the dimension of length. The 
characteristic frequency tu c of the fluctuations, that pro- 
vide the main contribution to the Casimir force, becomes 
smaller, 



if 



Id c 
o P \ - -C w p , uj c <C -, 
V a a 



d <C S = c/u)p, 5 <C a. 



(2) 



(3) 



This means that the high-temperature regime for the 
Casimir attraction of a film occurs at lower tempera- 
tures: at T > T c — huj c /k cc d 1 / 2 . In addition, under 
conditions ([3|), the surface impedance of a metal film is 
not small. Therefore, the Casimir force for thin films be- 
comes smaller than for bulk materials (see results of re- 
cent experiment [13| with thin films) , and the relative role 
of the temperature effects in the Casimir force becomes 
stronger. Thus, as we demonstrate below, the anomalous 
temperature dependence of the Casimir force can be ob- 
served, in principle, for thin metal films. The successful 
implementation of this experiment could put an end to 
the longstanding discussion on the role of the electron 
relaxation in the Casimir effect. 

Model. — The general formula for the Casimir interac- 
tion force between dielectric slabs with arbitrary dielec- 
tric constants e was orig inally derived by Lifshitz [3] 



(see, also, Refs. [15|, llq ). The Casimir force is pre- 



sented in this formula as a functional defined on the set 
of functions e(iu) n ) of a discrete variable uj n = 2imkT 
(n = 0, 1, 2, ... ). For the dielectric permittivity of the 
metal film, we choose the Drude dispersive model Eq. fl} 
which takes into account the temperature dependence of 
the relaxation frequency v caused by the scattering of 
electrons by phonons. We use the relation, 

v{T) = u + Up h (T/Q), 



v P h(x) = Av ph (l)x 5 



1/x 



y 5 dy 



(ev - 1) (l-e-J')' 



(4) 



based on the Gruneisen formula for the t emp erature de- 
pendence of the resistivity (see, e.g., Ref. \Ya) . Here v is 
the residual relaxation frequency caused by the electron 
scattering on crystal defects, Q is the Debye temperature, 
v p h(T/<d) is the relaxation frequency due to the electron- 
phonon scattering. The value ^ p h(l) depends on the 
Fermi velocity of electrons, the strength of the electron- 
phonon interaction, etc. This u p h(l) can be be obtained 
by measuring the resistivity at the Debye temperature. 

The constant A is ( y 5 dy/ (e y - 1) (1 - e~ v )j « 3. 
For simplicity, we do not take into account the surface 
scattering of electrons in the explicit form (Pfl) because it 
only changes the value of vq (see, e.g., Ref. Il8l). 

We consider first the Casimir effect for an ideal bulk 
conductor and a thin metal film of thickness d, separated 
by a distance a. Then, using the "Proximity Force Theo- 
rem" fl9| , we derive the expressions for the Casimir force 
between a metal film and an ideal metal sphere (lens). 
The geometry of problem is shown in Fig. 1. 





FIG. 1: (Color online) Geometry of the problem, (a) The 
Casimir attraction of a thin metal film (in red) to an ideal 
plane bulk metal, (b) The same film now interacting with a 
metal sphere. 

Casimir force. — The asymptotic equation for the 
Casimir attraction of a thin metal film to an ideal bulk 
plane metal was derived in Ref. [T3 . The force / per unit 
area can be written in the form, 



/ = - 



BkT 
8na 3 



dxx 3 e x I(x), 



where B = (htu c /4TrkT) 2 , u> c = u p (d/ 'a) 1 / 2 , 



OO 



1 



n(n + C) + B^(x) 



(5) 



(6) 



*(x) = x(l - e~ x ), C = hv(T)/2irkT, the prime over 
the sum symbol indicates that the term with n = is 
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taken with half the weight. Equation §5§ is valid if con- 
ditions ([3]) are fulfilled. In this case, one can neglect the 
relativistic retarding effect and pass to the limit c — > oo. 

Using the Abel-Plana formula for summing series, we 
can rewrite Eq. ([5]) in the form of a sum of two terms 
that correspond to two sources for the temperature de- 
pendence of the Casimir force, 



fu = ~ 



ftWr 



327r 2 a 3 



/ fv~^~ /rad 7 



dx x 3 t 



(It 



o t(t + rj) + 



fr&d 



hv{T) 
' 8w 2 a 3 j ( 
tdt 



dx : 



o (eT* -l){[t 2 -y{x)} 2 + V 2 t 2 }' 



(7) 



, (8) 



(9) 



where rj = 2v(T) /uj c , 7 = huj c /2kT. 

The first term in Eq. fTJ is provided by the quan- 
tum fluctuations of the electromagnetic field. It depends 
on temperature only via v(T). Obviously, the modulus 
of this term decreases when increasing the temperature. 
The asymptotics of /„ for low and high values of the pa- 
rameter 77 are: 



u = 

fu = - 



1 — «2 



647T0 3 

3hu> 2 



167r 2 a 3 i/(T) 



In 



U3> 



v<€.cj Ci (10) 
13 , v > w c , (11) 



where h w 3.51214, i 2 = 4C(3)/7rii « 0.43578, i 3 « 
0.59272, and is the zeta function. 

The second term in Eq. ([7]) is caused by the thermal 
fluctuations of the electromagnetic field. Its modulus in- 
creases when increasing the temperature. This term has 
different asymptotics in different temperature intervals: 



low-T 
rad 



at low temperatures and 



i/(T)(fcT) 2 fhuj c \ , _ _ , , 

(12) 



high-T 
rad 



C(3) kT 



kT ^> min (hu> c , hu 2 /v) (13) 



at high temperatures. In the case v <C cj c , there exist the 
intermediate asymptotics, 



intcrmcd — T 
rad 



C(3) (kT) 3 
2tt a 3 (huj c ) 2 



(14) 

Using Eqs. (|TU ]) -([l4" |) and the Proximity Force Theo- 
rem [l9(, one can easily derive the analogue asymptotics 
for the Casimir force F, 



F = 2-kR 



da'f(a'), 



(15) 



between an ideal metallic sphere of radius R and a thin 
metal film: 

F v (v < cj c ) = - 1 - — , (16) 



F^(zy>w c ) 



80a 2 



frw 2 i? 



4 w c 



87ra 2 i/(T) 

For the low-temperature interval in Eq. (Ti~2"j) 



i3 + c 




low-T 
rad 



TtRv(T){kTf 
12a 2 hu 2 
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for the high-temperature interval in Eq. 



high-T C(3) RkT 



ad 



(17) 



(18) 



(19) 



and for the intermediate temperature interval in Eq. (TT 



F 



intcrmcd — T 
rad 



-C(3) 



R(kT) 3 
a 2 (huj c ) 



(20) 



The above results show that the contribution to the 
Casimir force from quantum fluctuations always de- 
creases when increasing the temperature. At the same 
time, the term related to the thermal radiation always 
increases with temperature. As a result of this competi- 
tion, an interesting non-monotonous temperature depen- 
dence of the total Casimir force can be observed. Figure 
2 shows the dependence of the Casimir force between 
a metal film and: (a) an ideal metal semi-space, and 
(b) a metal sphere (lens), on temperature, calculated for 
different values of the parameter f pn (l). The Casimir 
force is normalized to the force -Fbuik between ideal bulk 
metals for the same separation a. The force -Fbulk is, 



F 



plane 



bulk 



7r he S /240a , for bulk samples with planar 



boundaries of the area S, and F^ le = 7r 3 /ici?/360a 3 for 
the attraction between a metal semi-space and a metal 
sphere with radius R. Note that wc use the relations 
for -Fbuik in the retardation limit because c/uj p a < 1. 
For the same separations a, the Casimir force for a thin 
film is calculated for the non- retardation regime, since 
c/uj c a = c/uj p (da) 1 / 2 3> 1. 

The decreasing portion of the F(T) dependence corre- 
sponds to a decrease of the Casimir contribution to the 
entropy when increasing T. This decrease is connected to 
the enhancement of the electron scattering on phonons 
and is much weaker than the increase of the phonon con- 
tribution to the entropy. Thus, the total entropy cer- 
tainly increases when increasing the temperature. 

In conclusion, we predict an unusual non-monotonous 
temperature dependence of the Casimir attraction force 
between a thin metal film and a bulk plane ideal metal 
or a metal sphere (lens). Usually, for bulk samples, the 
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FIG. 2: The temperature dependence of the Casimir force of 
attraction of a metal film to an ideal plane bulk metal (a) 
and to a metal sphere (b) for different separations a. The 
left vertical axes show F (T)/i<buik and the right vertical axes 
show 1 — F(T) / F(0). Curves I, 2, 3, 4, 5 correspond to u p h • 
(lCT 13 s) = 0.7, 1.05, 1.4, 1.75, 2.1, respectively. The values 
of the other parameters are: u> p — 2 ■ 10 15 s _1 , vo = 10 11 s -1 , 
9 = 100 K. 

Casimir force increases slowly with temperature. Here 
we predict a noticeable decrease of the force with an in- 
crease of T for metal films. The experimental observation 
of this unusual temperature dependence of the Casimir 
force can put an end to the long-standing dispute on the 
role of the electron relaxation in the Casimir effect. 

We acknowledge partial support from the NSA, LPS, 
ARO, NSF grant No. EIA-0130383, JSPS-RFBR 06-02- 
91200, MEXT Grant-in-Aid No. 18740224, the EPSRC 
via No. EP/D072581/1, EP/F005482/1, ESF AQDJJ 
network programme, and the JSPS CTC Program. 



[1] P.W. Milonni, The Quantum Vacuum (Academic press, 

San Diego, 1994). 
[2] K.A. Milton, The Casimir Effect (World Scientific, Sin- 

gapure, 2001). 

[3] M. Kardar and R. Golestanian, Rev. Mod. Phys. 71, 1233 
(1999). 



[4] F. Capasso, J.N. Munday, D. Iannuzi, and H.B. Chan, 
IEEE J. Sel. Top. Quant. El. 13, 400 (2007). 

[5] C. Cattuto, R. Brito, U. M. Marconi, F. Nori, and 
R. Soto, Phys. Rev. Lett. 96, 178001 (2006). 

[6] S.K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997); U. Mo- 
hideen and A. Roy, Phys. Rev. Lett. 81, 4549 (1998); F. 
Chen et al, Phys. Rev. A 69, 022117 (2004); R.S. Decca 
et al, Ann. Phys. (NY) 318, 37 (2005); D. Iannuzzi et 
al, Solid State Commun. 135, 618 (2005); J.N. Munday 
et al, Phys. Rev. A 71, 042102 (2005); J.N. Munday et 
al, New J. Phys. 8, 244 (2006); D. Iannuzzi et al, J. 
Phys. A: Math. Gen. 39, 6445 (2006); J.N. Munday and 
F. Capasso, Phys. Rev. A 75, 060102(R) (2007). 

[7] H.B. Chan et al, Science 291, 1941 (2001); H.B. Chan 
et al, Phys. Rev. Lett. 87, 211801 (2001). 

[8] M. Bordag et al, Phys. Rev. Lett. 85, 503 (2000); 
M. Bordag et al, Phys. Rev. Lett. 87, 259102 (2001); 
V.B. Bezerra, G.L. Klimchitskaya, and V.M. Mostepa- 
nenko, Phys. Rev. A 66, 062112 (2002); F. Chen et al, 
Phys. Rev. Lett. 90, 160404 (2003). 

[9] M. Bostrom and Bo E. Sernelius, Phys. Rev. Lett. 84, 
4757 (2000). 

[10] Bo E. Sernelius and M. Bostrom, Phys. Rev. Lett. 87, 
259101 (2001); I. Brevik, J.B. Aarseth, and J.S. Hoye, 
Int. J. Mod. Phys. A 17, 776 (2002); I. Brevik, 
J.B. Aarseth, and J.S. Hoye, Phys. Rev. E 66, 026119 

(2002) ; J.S. Hoye, I. Brevik, J.B. Aarseth, and K.A. Mil- 
ton, Phys. Rev. E 67, 056116 (2003); C. Genet, A. Lam- 
brecht, and S. Reynaud, Phys. Rev. A 67, 043811 

(2003) ; I. Brevik, J.B. Aarseth, J.S. Hoye, and K.A. Mil- 
ton, Phys. Rev. E 71, 056101 (2005); I. Brevik and 
J.B. Aarseth, J. Phys. A: Math. Gen. 39, 6589 (2006); 
I. Brevik, S.A. Ellingsen, and K.A. Milton, New J. Phys. 
8, 236 (2006); I. Brevik, S.A. Ellingsen, J.S. Hoye, and 
K.A. Milton, Preprint arXiv:0710.4882vl [quant-ph] 25 
Oct 2007. 

[11] E.M. Lifshitz and LP. Pitaevskii, Statistical Physics, 
Part 2 (Pergamon Press, Oxford, 1980), Sec. 81. 

[12] V.N. Dubrava and V.A. Yampol'skii, J. Phys. A: Math. 
Gen. 33, L243 (2000). 

[13] M. Lisanti et al, Proc. Nat. Acad. Sci. USA 102, 11989 
(2005). 

[14] E.M. Lifshitz, JETP 2, 73 (1956). 

[15] J. Schwinger, L.L. DeRaad, and K.A. Milton, Ann. Phys. 

(N.Y.) 115, 1 (1978). 
[16] J. Schwinger, Proc. Nat. Acad. Sci. USA 89, 4091 (1992). 
[17] H. Ibach and H. Liith, Solid-State Physics, (Springer, 

Berlin, 2003). 

[18] N.M. Makarov, A.V. Moroz, and V.A. Yampol'skii, Phys. 

Rev. B 52, 6087 (1995). 
[19] J. Blocki, J. Randrup, W.J. Swiatecki, and C.F. Tsang, 

Ann. Phys. (N.Y.) 105, 427 (1977). 



